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ABSTRACT 


Two  previous  papers^^  have  described  the  concept  and  some  of 
the  applications  of  Narrowband  Interferometry,  using  monostatic  and 
multistatic  radars  (at  small  bistatic  angles)  to  obtain  two  and  three 
dimensional  "images"  of  targets  in  terms  of  their  scattering  centers. 

The  method  involves  taking  finite  Fourier  transforms  of  measured  narrow- 
band  amplitude  and  phase  data  to  resolve  individual  scattering  centers 
in  "Doppler"  at  each  radar  site. 

The  present  paper  extends  these  methods  in  two  directions.  The 
class  of  scattering  centers  considered  is  widened  to  include  not  only 
isotropic  point  scatterers  as  before  but  also  point  dipoles,  as 
representatives  of  orientable  scattering  centers.  The  theory  is 
broadened  to  involve  the  Polarization  Scattering  Matrix,  so  that  the 
polarization  behavior  of  the  individual  resolved  scatterers  can  be 
determined  as  aids  to  their  recognition. 
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SECTION  I 


INTRODUCTION 


Two  previous  papers^  have  described  the  concept  and  some  of 

the  applications  of  Narrowband  Interferometry  using  monostatic  and  multi¬ 
static  radars  (at  small  bistatic  angles)  to  obtain  two  and  three  dimen¬ 
sional  "images"  of  targets  in  terns  of  their  scattering  centers.  The 
systems  therein  discussed  did  not  use,  nor  did  the  theory  refer  to, 
measurements  made  in  more  than  one  polarization.  This  was  sufficient 
for  studying  objects  composed  only  of  simple  isotropic  point  scattering 
centers.  But  real  objects  do  not  generally  contain  many  such  simple 
scatterers,  and  those  that  do  exist  are  of  low  cross  section  compared 
with  other  types  of  scattering  centers.  In  particular,  realistic  targets 
are  better  described  in  terms  of  a  model  which  allows  for  orientable 
"scatterers1',  since  such  things  as  edges  are  usually  pre-eminent  in  the 
scattering-center  description  of  a  real  object.  But,  in  the  case  of 
orientable  scattering  centers,  valuable  information  is  contained  in  the 
terms  of  the  complete  Polarization  Scattering  Matrix,  which  information 
will  be  lost  if  the  Scattering  Matrix  is  not  utilized. 

Now,  under  ARPA  sponsorship,  The  MITRE  Corporation  has  recently 
developed  and  installed  a  radar  system  capable  of  measuring  the  full 
scattering  matrix  on  a  pulse-by-pulse  basis;  this  system  has  been  inte¬ 
grated  into  the  MITRE-Lincoln  Laboratory  Tristatic  Radar  Interferometer. 
Therefore,  the  present  paper  will  be  devoted  to  considering  what  informa¬ 
tion  is  present  in  these  scattering  matrix  measurements  and  in  particular 
how  an  appropriate  Fourier  Transformation  of  the  terms  of  the  scattering 
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matrix  can  help  to  identify  the  polarization  properties  of  those  indi¬ 
vidual  scattering  centers  which  can  be  resolved  in  this  manner.  The 
theory  will  be  devoted  primarily  to  two  types  of  scattering  centers: 
isotropic  point  scatterers  as  simple  representatives  of  non-or ientab le 
scattering  centers,  and  point  dipole  scatterers  as  simple  representatives 
of  orientable  scattering  centers.  Although  these  idealized  scattering 
centers  are  in  some  respects  considerable  abstractions  from  the  more 
complicated  sorts  of  scatterers  found  on  real  targets,  they  have  the 
useful  property  and  great  advantage  of  illustrating  most  of  the  pertinent 
geometric  properties  of  the  more  complex  scatterers  without  the  latter's 
attendent  mathematical  complexities  which  obfuscate  without  enlightening. 

(3) 

As  a  matter  of  fact,  a  recent  study  of  the  "Doppler  Mapping"  of  an  Edge 
using  the  methods  of  the  present  paper  has  clearly  borne  out  the  contention 
that  the  point  dipole  is  a  good  approximation  to  the  more  realistic  edge 
scatterer . 

In  connection  with  the  term  "Doppler  Mapping','  given  in  quotation 
marks  in  the  previous  paragraph,  I  feel  that  scxne  comments  are  deserved. 

The  process  of  taking  finite  Fourier  Transforms  of  amplitude  and  phase 
data  in  the  time  domain  has  been  called  "Doppler  Processing"  or  "Doppler 
Mapping"  here  at  MITRE  for  a  long  time.  The  resulting  "Doppler  Maps" 
are  said  to  show  scattering  centers  resolved  "in  Doppler",  by  which  is 
meant  resolution  in  the  f  (frequency)  domain.  But  it  will  became  quite 
apparent,  upon  careful  examination  of  the  analysis  given  in  Section  HI. 
that  these  terms  are  misnomers.  Although  the  mathematical  forms  of  the 
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functions  obtained  are  identical  to  those  which  are  characteristic  of  the 
Doppler  Effect,  the  fact  is  that  the  terms  due  to  the  true  Doppler  effect 
are  negligible  compared  to  other  terms,  and  drop  out  of  all  the  equations 
at  an  early  stage.  The  effect  responsible  for  the  resolution  of  scatter¬ 
ing  centers  in  frequency  is  due  to  periodic  changes  in  the  phase  differ¬ 
ences  between  different  scattering  centers  as  these  centers  move  with 
respect  to  one  another.  Nevertheless,  the  erroneous  terms  "Doppler  Map", 
etc.,  have  become  standard  through  long  usage  here  at  MITRE,  and  they 
will  continue  to  be  used  in  this  paper,  though  with  quotation  marks  wherever 
I  could  remember  to  put  them  in.  Whenever  quotation  marks  do  not  appear 
on  these  terms,  the  reader  is  asked  please  to  imagine  them. 
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SECTION  II 


THE  POLARIZATION  SCATTERING  MATRIX 


The  received  signal,  in  an  arbitrary  polarization,  is  related 
to  the  transmitted  signal,  in  an  arbitrary  polarization,  through  the 
Polarization  Scattering  Matrix.  The  form  of  the  matrix  depends  on  the 
modes  of  polarization  being  used  and  on  the  representation  employed  to 
describe  these  modes.  We  will  be  concerned  here  with  two  different  modes, 
linear  and  circular.  As  for  the  representations,  a  glance  at  the  litera¬ 
ture  will  show  that  there  is  as  yet  no  generally  accepted  convention, 
but  that  the  situation  is  in  a  state  of  great  confusion.  After  reading 
the  interesting  discussions  of  the  Scattering  Matrix  to  be  found  in 
References  4,  5,  6,  and  7,  and  giving  great  weight  to  the  conventions 
currently  being  employed  by  those  directly  involved  in  the  ARPA  Scattering 
Matrix  Radar  Program  here  at  MITRE,  and  with  due  regard  to  the  additional 
requirement  that  the  form  of  the  scattering  matrix  for  a  bistatic  situation 
reduce  to  that  for  the  monostatic  case  when  the  bistatic  angle  goes  to 
zero,  I  have  chosen  to  define  things  in  what  I  hope  is  a  reasonable 
manner.  Any  differences  between  the  present  representation  and  those 
found  in  the  above  references  (aside  from  those  caused  merely  by  the 
choice  of  which  polarizations  are  denoted  1  and  2)  will  be  in  the  signs 
of  the  phase  factors  and  should  not  cause  any  significant  problems. 
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1.  Linear  Polarization 


Consider  a  bistatic  situation  in  which  a  linearly  polarized 

A 

wave  is  transmitted  in  direction  k,  strikes  a  target,  0,  in  the  far- 

A 

field  and  is  reflected  toward  a  receiver  in  direction  -k'  relative  to 
the  target;  see  Figure  1. 


Figure  1.  The  Bistatic  Situation 


The  bistatic  angle,  p,  is  defined  by 

A  /V 

cos  (3  =  k  •  k'  , 

and  unit  vectors  n  and  n'  are  normals  to  the  earth's  surface  at  the 
transmitter  and  receiver  respectively. 

A  plane  normal  to  k  (or  k')  is  called  the  Transmitter  (or  Receiver) 
Polarization  plane.  In  the  transmitter  polarization  plane  we  may  define 
a  unit  vector,  ex  ,  by  the  equation 


=  n  x  k 

|n  X  k| 


(2.1a) 
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Ihis  vector  is  clearly  'horizontal"  according  to  the  usual  meaning  of  the 
term,  and  a  datum  line  parallel  to  this  vector  in  the  transmitting  antenna 
plane  may  be  maintained  always  horizontal  however  the  antenna  rotates,  thus 
providing  a  fixed  datum  for  "horizontal  polarization"  in  transmission. 

We  can  now  define  another  unit  vector,  ea  ,  in  the  transmitter  polarization 
plane,  such  that  (ej ,  ea  ,  k)  form  a  right-handed  ortho-normal  triad: 


ea  =  k  X  ex  (2.1b) 

A  datum  line  parallel  to  this  vector  in  the  transmitting  antenna  plane 
is  always  normal  to  the  "horizontal  polarization"  datum  line  and  thus 
provides  a  fixed  datum  for  what  is  usually  called  "vertical  polarization" 
in  transmission,  even  though  it  is  not  usually  a  "vertical"  line  in  the 
common  meaning  of  that  term. 

Similarly,  to  represent  the  polarization  vector  of  the  received 
signal  we  choose  basis  vectors  e^ '  and  ea  1  in  the  receiver  polarization 
plane: 


n'  X  k' 

| n *  X  k'| 


(2.2a) 


ea  =  k'  X  ej[  (2.2b) 

Once  again,  the  respective  datum  lines  in  the  receiving  antenna  plane 
provide  fixed  lines  for  defining  "horizontal"  and  "vertical"  polarization 
of  the  received  signal.  Moreover,  the  definitions  above  guarantee  that, 
since  k  k  and  n1  -*  n  as  p  -+  0,  these  two  sets  of  bases  vectors  (and 

lines)  have  the  desired  property  of  becoming  identical  as  (3  -♦  0. 
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Now  we  must  exert  just  a  bit  of  care  to  define  the  Scattering 
Matrix  in  a  way  that  is  valid  for  both  the  bistatic  and  monostatic  cases. 
Suppose  the  transmitting  antenna  sends  out  a  unit  signal  with  an  arbitrary 
linear  polarization,  represented  naturally  in  the  unprimed  basis  as 


ET  =  ET  =  a1e1  +  aaea 


(2.3) 


where  (a^2  +  (aa)2  =  1.  After  striking  the  target  and  starting  toward 
the  receiver,  the  returning  signal,  now  represented  naturally  in  the  primed 
basis,  may  be  written  as 


+  C 


a  e2 


(2.4) 


We  can  now  define  the  scattering  Matrix  S  by  requiring  that 


eR  =  S  Et  .  (2.5) 

This  definition  differs  from  the  usual  definition  of  the  scattering  matrix 
in  that  it  connects  two  polarization  vectors  expressed  in  different  co¬ 
ordinate  systems;  however,  it  reduces  to  the  usual  definition  when  (3-*0. 

In  order  to  measure  the  elements  of  S,  we  assume  that  the  receiving 
antenna  is  linearly  polarized  in  some  direction  E^  which  is  represented 
in  the  primed  basis  as 
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1.  Then  the  signal  which  is  actually  received  is 


where  (a^)2  +  (a^)2  = 


(2.7) 


This  tells  us  that  in  order  to  measure  element  of  the  scattering  matrix 
we  must  transmit  a  unit  signal  with  polarization  e^  and  measure  the  received 
signal  with  polarization  e|  ,  whereupon  the  value  of  the  received  signal 
is  equivalent  to  the  element  . 


2.  Circular  Polarization 

It  is  with  circular  polarization  that  the  confusion  of  represent¬ 
ation  referred  to  above  enters  to  becloud  the  picture.  We  will  have  to 
make  several  definite  choices,  which  will  be  indicated  as  we  go  along. 

One  of  these  occurs  in  the  definition  of  "Right  Circular"  and  "Left 
Circular"  polarization.  We  will  not  sow  confusion  by  describing  the 

conventions  which  will  NOT  be  used,  but  rather  will  describe  only  the 

_  ( 8 ) 

standard  IEEE  convention  which  WILL  be  used.  This  convention  defines 

the  "handedness'1  of  the  circular  polarization  in  term#  of  the  direction 

of  propagation  of  the  signal:  A  right-circularly  polarized  wave  is  one 

for  which  the  E-vector  rotates  in  tha  direction  of  a  right-handed  screw 

as  it  advances,  and  conversely  for  a  left-circular ly  polarized  wave.  To 


put  it  more  graphically,  point  the  thumb  in  the  direction  of  propagation 

so  that  the  fingers  curl  in  the  direction  of  rotation  of  the  E-vector; 

then  the  "handedness"  of  the  circular  polarization  is  defined  by  the  physical 

hand  which  must  be  used  to  accomplish  this.  , 

Let  us  now  define  a  basis  for  representing  circular  polarization. 

The  unit  vectors  e1}  ez,  £  as  defined  above  form  a  right-handed  orthonormal 

triad,  and  the  same  goes  for  their  primed  counterparts.  Consider  a  plane 

* 

wave  advancing  along  the  positive  k  direction  with  polarization  vector 
given  by  ex  +  ie2  .  Of  course  what  we  really  mean  by  this  is  a  real  wave 
of  the  form 

E  =  Re  {(e1+iea)e1^’r”UJt^  }  =  ex cos (k* r-ujt)  -  e2 sin(k* r-uut)  . 

Examining  the  angle  0  between  the  polarization  vector  and  the  ex  axis 
we  see  that 


tan  0 


sin(k»r-uut) 

cos(k*  r-uut) 


tan(uut-k'  r) 


whence 


0(r,t)  *  LUt-k * r  ; 

this  shows  that  as  a  function  of  t  the  polarization  vector  is  advancing 
and  turning  in  the  same  way  as  a  right-handed  screw. 
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Next  let  us  consider  a  plane  wave  advancing  along  the  (-k)  direction 
with  the  same  polarization  vector  e1  +  iea  .  In  this  case  what  is  meant 
is  a  real  wave  of  the  form 


E  H  Re  {(e1+ie2)e1^k’r+UJt^  }  =  ex  cos (k*  rfO)t)  -  e2  sin(k- r+uut)  , 


whereupon  the  angle  9  between  the  polarization  vector  and  the  el  axis 
becomes 


6(r  ,t) 


-uut-k*  r 


Since  the  wave  is  propagating  in  the  direction  (-k)  ,  this  shows  that  as 
a  function  of  t  the  polarization  vector  is  again  advancing  and  turning 
in  the  same  way  as  a  right-handed  screw. 

Thus  we  have  seen  that  the  complex  polarization  vector  ex+ie2 
always  represents  a  r ight-c ircularly  polarized  wave,  both  in  transmission 
and  reflection.  Similarly,  we  could  show  that  the  complex  polarization 
vector  ej^-i^  always  represents  a  lef t-circularly  polarized  wave.  There¬ 
fore  we  define,  as  basis  vectors  for  representing  circular  polarization, 

A  A 

the  complex  vectors  and  Cg  which  in  the  linear  basis  have  the 
representat ions 


J2 


Cl  ■ 


ei 


+  iea  = 


R.H.C. 


(2.8a) 


L.H.C. 


(2.8b) 
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A  similiar  basis  is  defined  for  the  primed  coordinate  system. 


A  word  about  normalization  may  be  in  order  here.  Because  they 
are  complex,  the  vectors  cannot  be  normalized  in  the  usual  dot  product 

sense.  Instead  we  must  use  the  Hermit ian  scalar  product  of  the  form 

-*  -•  —  —  * 

<  U,V  >  =  U*V  , 

whereupon  we  see  that 


<  C 


i’ 


V 


tj 


But  actually  it  will  not  be  necessary  to  use  this  form  of  normalization 

in  the  analyses  to  follow;  in  use  the  always  appear  multiplying  a 

"phase  factor"  of  the  form  .  .r*  -•  ,  . 

ei(k*r  +  u)t) 

with  the  Real  Part  of  the  product  being  understood  to  be  taken.  When 
this  has  been  done,  the  resulting  vectors  will  be  found  to  be  properly 
normalized  in  the  ordinary  sense,  so  no  more  need  be  thought  about  it. 

Let  us  now  proceed  to  transform  the  scattering  matrix  from  its 
representation  in  the  linear  bases,  S,  to  its  representation  in  the 
circular  bases,  which  we  will  denote  by  C  .  Let  A  be  the  matrix  that 
transforms  a  vector  from  the  circular  basis  to  the  linear  one;  using 
eqns.  (2.8) it  is  easily  seen  that 

A  -  >  p  ll 

v2  Li  -ij  (2.9) 
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Now  we  rewrite  eq.  (2.7)  in  the  form 


V  =  ER(e  1 )  S  ET(e),  (2.10) 

where  the  argument  e1  or  e  refers  to  the  primed  or  unprimed  linear 
basis  and  the  tilda  (~)  indicates  the  transpose  of  a  vector  or  matrix  . 
In  an  obvious  notation, 


ER(e')  =  A  ER(c ')  and  E^(e)  =  AET(c)  , 


so  eq.  (2.10)  becomes 


V  *  Er(c')A  SAEt(c)  .  (2.11) 

But  the  scattering  matrix  C  is  defined,  analogously  to  the  definition 
of  S ,  by  the  identity 


V  =  Er(c ' ) *  C  Et  (c)  ,  (2.12) 

whereupon  we  see  that 

C  =  A  s  A  (2.13) 
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A  little  matrix  multiplication  now  yields  the  elements  of  C: 


2C11  - 

(sn- 

s22)  + 

1<S12+ 

S21 

2C12  " 

(sn+ 

s22) 

‘<s12- 

S21 

2C21 

<su+ 

s22>  ♦ 

CN| 

f—i 

C/5 

V-/ 

•H 

S21 

2C22  " 

(sn- 

S22)  - 

i(S12+ 

S21 

(2.14) 


We  should  also  consider  the  effect  of  Faraday  rotation,  which 
takes  a  simple  form  in  the  circular  representation.  Faraday  rotation 
is  due  to  an  anisotropy  of  the  ionosphere  which  causes  the  propagation 
constant  for  an  electromagnetic  wave  to  be  slightly  different  according 

(4) 

to  the  sense  of  circular  polarization  of  the  wave  traversing  the  medium. 

The  result  is  that  one  component  gains  in  phase  proportional  to  the  distance 
traversed,  and  the  other  component  loses  an  equal  amount  of  phase.  The 
discussion  of  this  subject  in  Reference  4  is  very  clear  and  need  not 
be  repeated  here;  the  net  result  is  that  the  scattering  matrix  C  has 
additional  phase  factors  in  the  off-diagonal  terms.  Viz: 


'Cll 

C12 

C  e'12’*' 

C21e 

C22 

+i2i|f 


(2.15) 


where  ,  the  amount  of  one-way  Faraday  rotation,  is  proportional  to  the 

distance  between  the  radar  and  the  target.  Although  it  varies  with  time 
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of  day,  the  sunspot  cycle,  etc.,  the  value  of  \|f  for  a  traversal  of  the 
entire  ionosphere  has  not  exceeded  30°  in  measurements  made  at  MITRE 
it  is  generally  a  rather  slowly  varying  function  of  distance  and  time. 


M 


SECTION  III 


THE  SCATTERING  CENTER  MODEL 

In  this  paper  we  adopt  the  model  whereby  the  scattering  from  an 

object  is  presumed  to  be  explained  in  terms  of  a  relatively  small  number 

of  isolated  scattering  centers  of  various  sorts;  the  description  of  the 

object  then  consists  of  determining  the  locations  and  properties  of  its 

constituant  scattering  centers.  Some  of  these  scattering  centers  (e.g. 

point  scatterers)  are  assumed  to  be  rigidly  attached  to  the  center  of 

rotation  of  the  object,  so  that  a  study  of  the  motions  of  these  scatterers 

enables  one  to  determine  the  gross  body  motion  of  the  object  itself. 

Others  (e.g.  edges)  are  only  in  part  rigidly  connected  to  the  object 

center,  being  capable  of  "sliding"  along  a  space  curve  as  the  object 

rotates;  such  scatterers  are  rather  more  difficult  to  deal  with,  and  pro- 

(3) 

vide  less  direct  information  about  the  gross  body  motion.  Yet  another 
type  of  scattering  "center",  the  specular,  is  slippery  in  two  dimensions, 
being  capable  of  "sliding"  along  a  surface  in  space  as  the  object  rotates; 
this  type  partakes  sb  little  of  the  characteristic  of  localizability  that 
perhaps  it  ought  to  be  excluded  from  the  class  of  "scattering  centers" 
altogether. 

The  scattering  center  model  of  an  object  has  been  in  use  for  quite 

some  time,  and  it  has  been  employed  in  numerous  problems.  It  appears  to 

be  the  only  convenient  way  to  deal  with  the  subject  of  the  present  paper, 

and  it  will  be  used  herein.  But,  though  this  is  not  the  time  or  place 

for  a  critical  review  of  the  concept,  a  serious  caveat  is  definitely  called- 

for.  The  model  must  be  understood  in  its  context  as  a  high  frequency 
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approximation,  that  is,  that  the  object  being  studied  and  in  fact  the 
individual  scattering  centers  themselves  are  large  compared  with  the 
wave-length.  In  this  context  it  makes  little  sense  to  speak  of  small 
isotropic  "point"  scatterers:  the  nearest  thing  to  a  true  rigidly-attached 
point  scatterer  is  a  tip,  which  is  approximately  isotropic  but  is  very 
low  in  reflecting  power;  on  the  other  hand,  the  only  sizeable  truely  iso¬ 
tropic  scattering  center  is  a  specular,  which  also  most  resembles  a  "point" 
but  is  far  from  being  rigidly  attached.  All  the  other  types  of  scattering 
centers  are  also  more  or  less  idealizations  of  the  true  situation.  The 
most  realistic  is  the  edge,  which  also  probably  makes  up  the  class  of 
scatterers  most  prevalent  in  objects  other  than  smooth  ellipsoids.' 

What  is  one  to  infer  from  all  this?  We  must  deal  with  a  mathemati¬ 
cally  tractable  model,  or  else  the  complexities  of  the  solution  will  obscure 
the  simple  physical  realities  which  are  being  presented.  So  we  will  use 
isotropic  point  scatterers,  and  we  will  represent  orientable  scatterers 
(such  as  edges)  primarily  by  the  fiction  of  point  dipoles,  and  we  will 
touch  on  real  edges  only  a  slight  bit  (though  in  fact  their  properties 
seem  to  be  well  represented  in  terms  of  point  dipoles  that  can  "slip") , 
but  we  will  also  keep  in  mind  that  these  are  models ,  whose  only  true  test 
of  validity  is  this:  does  it  help  to  explain  the  experimental  observations? 
To  the  extent  that  it  does,  the  model  (however  unrealistic)  is  justified; 
to  the  extent  that  it  does  not,  we  must  seek  for  a  more  reasonable  mode1. 

With  this  understanding,  then,  let  us  describe  the  scattering  center 
model  and  its  relationship  with  the  radar  system,  and  consider  carefully 
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what  happens  as  a  transmitted  radar  signal  strikes  the  target  and  returns 

to  the  receiving  antenna.  We  assume  that  the  target  is  composed  of  a  finite 

(and  relatively  small)  number  of  scattering  centers  of  various  kinds,  all 

connected  in  some  way  to  a  center  of  mass  0  which  is  being  translated 
— ♦ 

with  velocity  V(t)  .  The  target  is  supposed  to  be  rotating  about  0 
as  a  rigid  body  with  angular  velocity  U)(t)  .  Some  of  the  scattering  centers 
(such  as  the  idealized  isotropic  point  scatterers  and  point  dipoles)  are 
assumed  to  be  rigidly  connected  to  the  center  of  mass  and  therefore  par¬ 
take  fully  of  the  rigid  rotation  uu  .  Others  (such  as  edges  and  speculars) 
must  be  visualized  as  slipping  along  a  space  curve  or  a  surface  and  hence 
partake  only  in  part  of  the  rigid-body  rotation.  But  all  of  them,' whether 

slippery  or  not,  may  be  assumed  to  be  located  at  a  particular  point  at 

— ♦ 

each  instant  of  time,  given  by  r^(t)  with  respect  to  0. 

This  target  is  being  tracked  by  a  radar  interferometer  consisting 
of  a  narrowband  transmitter-receiver  located  at  A  and  remote  receivers 
located  at  B  and  C;  the  angle  subtended  at  the  target  by  the  interfero¬ 
meter  is  assumed  to  be  small.  It  is  further  assumed  that  the  transmitter 
is  capable  of  transmitting  and  each  receiver  is  capable  of  receiving 
signals  of  two  orthogonal  polarizations,  and  each  receiver  can  measure 
both  amplitude  and  phase  in  each  polarization  channel.  The  radar  operates 
at  frequency  f^  and  transmits  pulses  of  length  A  seconds,  where  for 
concreteness  in  visualization  we  may  use  values  like  those  of  the  present 
MITRE-Lincoln  Laboratory  Interferometer,  i.e. , 
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f  =“  10^  MHz  and  A  =**  1  ms  . 
o 

Referring  to  Figure  2,  the  line-of-sight  vectors  from  radars  A,  B, 

— # 

,  R^.  For  all  reasonable 
satellite  objects  and  passes  the  following  set  of  conditions  will  easily 
obtain: 

1.  The  target  is  very  small  compared  with  its  distance  from  the 
radars;  in  fact 


C,  to  target  center  0  are  respectively  R^,  R^ 


max  {  I  r  -  r 
u  1  n  m 


min 


{R} 


I)  -4 
—  <  10 


2.  The  target  is  very  small  compared  with  the  length  of  the  radar 
pulse  in  space;  in  fact 


max  { I r  -  r  I } 
_  1  n  m1 

c  A 


<  10 


3. 

is  small; 


The  bistatic  angle  at  the  target  between  any  pair  of  radars 
in  fact 


p  ~  D/r  £  io'1 


radians . 
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Scattered 


v\ 


Figure  2.  Geometry  for  a  Two-Site  Interferometer 

(not  to  scale) . 
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* 


Now  assume  that  at  time  t^  antenna  A  transmits  a  pulse,  of 

A 

polarization  E^,  toward  the  target.  Conditions  1  and  2  above  allow 
us  to  make  the  far  field  approximation  whereby  the  transmitted  pulse  as 
it  approaches  scattering  center  n  may  be  considered  to  be  of  the  form, 


~  E  exp(2TTi£f!  [<ItA  +  •  r  )  -  c(c-Co>  ]}  ,  (3.1) 

ra 


where 


f 

n 


1  - 

-  Ra- 
c  A 


(3.2) 


represents  the  Doppler-shifted  frequency  of  the  incoming  wave  relative  to 

the  scatterer;  here  uu^  is  the  effective  angular  velocity  of  the  scatterer 

taking  into  account  any  slipperiness:  for  rigidly  attached  isotropic  points 

and  point  dipoles  =  uu.  In  writing  eq.  (3.1)  we  have  ignored  a  term  of 

order  r  /R  in  the  denominator  and  also  ignored  a  difference  in  angle 
n  a 

of  the  order  of  r^/R^  in  the  polarization  vector  at  the  location  of 

A 

scatterer  n  relative  to  E^,(incident  at  0);  both  of  these  approximations 
are  of  course  well  justified  by  condition  1.  The  use  of  a  monochromatic 
spherical  wave  for  the  incident  pulse  near  the  target  is  justified  by 
condition  2. 
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Notice  that  the  incident  wave  front  reaches  each  scatterer  n  at 

x 

location  xn  =  •  r^  at  time  t^  =  tQ  +  ~  ,  whereupon  eq.  (3.1) 

tells  us  that 

E  ~  Et  e‘®°/RA 

at  each  scattering  center  at  the  instant  it  begins  to  be  excited  by  the 
incident  pulse.  The  effect  of  the  scattering  center  on  the  incident 
pulse  will  be  symbolically  represented  by  the  operator  a  ,  so  that  the 
wave  which  the  scattering  center  begins  to  re-radiate  has  the  form 
aETe  ^°/ra  at  the  instant  when  re-radiation  begins.  The  form  of <the 
operator  a  will  of  course  depend  on  the  nature  of  the  scattering  center 
involved.  In  re-radiation  two  cases  must  now  be  distinguished. 

CASE  I.  RE-RADIATION  BACK  IN  THE  DIRECTION  OF  TRANSMITTER  A. 

The  wave  re-radiated  by  scattering  center  n  goes  in  the  direction 

A 

-R^  .  The  Doppler  shift  causes  an  additional  change  of  frequency  such  that 
the  new  frequency  is 


f* 


A  9  A 

R.  -  -  R.  • 
A  c  A 


UU  X 

n 


(3.3) 


Now  as  it  approaches  the  receiver  at  A  the  re-radiated  signal  wave  front 
has  the  form 
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r?A 

E  ~ 
n 


a  E„ 


exp  {2ni  fA 

r ,2 D 


t>  V 


rn+  c(t-tn)]) 


Here  t  is  the  nominal  time  of  arrival  of  the  re-radiated  wave  front  from 
the  center  of  the  target  at  0,  given  by  t  =  t  +  2R  /c;  we  use  this  nominal 

O  A 

time  so  that  we  can  later  sum  the  contributions  from  all  the  scattering 
centers  at  the  same  instant  of  reception.  With  this  value  of  t  the 
received  signal  from  scatterer  n  appears  as 


n 


CT  Et  e  exp  [-4rri  n  (R^  * 
A  c 


(3.4) 


Using  eq.  (3.3)  the  bracketed  phase  factor  in  eq.  (3.4)  can  be 
written  as 


0 


2 

c 


V-R.+ 

A 


2 

c 


R.X 

A 


But  for  reasonable  objects  and  passes  r  <  10m.,  uu  <  1  rad/sec,  and 

n  n 

V  ~  104m/sec,  whereupon 


and 


<  10’4 
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This  means  that  even  when  multiplying  the  factor  —  these  two 

o 

terms  contribute  no  more  than  about  3  to  the  phase  , and  we  shall  neglect 
them  in  comparison  with  the  first  term.  Then  we  may  write  the  received 
signal  from  scattering  center  n  as 


— *A  e  i(^0  /v 

En~—  aETe*P 

A 

If  the  receiving  antenna  is  polarized  in  direction  E  [where  E  •  R  =  0] 

R  R  A 

then  the  received  signal  from  scatterer  n  has  the  form 


ef  v  ^ 


e^o  iA 

W~  ER 

A 


a  ETexp 


r  ^TT  ~ 

r  iRi 

O 


(3.5) 


There  is  nothing  unexpected  about  this  result,  but  in  the  first  place  it 
demonstrates  that  there  is  no  direct  effect  from  the  Doppler  shift  on  the 
phase  of  the  signal  received,  and  in  the  second  place  it  shows  how  to 
proceed  in  the  less  obvious  bistatic  case  to  be  considered  next. 


CASE  II.  RE-RADIATION  IN  THE  DIRECTION  OF  REMOTE  RECEIVER  B. 

The  wave  re-radiated  by  scattering  center  n  goes  in  the  direction 
-R^  .  The  Doppler  shift  causes  an  additional  change  of  frequency  such 
that  the  new  frequency  is 


fn  “  4 [1  -  i  <V  V  -  i  <V  V  •  V  ^ 
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Now  we  define 


and 


ra  +  bb 
l*A  +  *bI 


cos  8  *  Ra-  Sj  , 


whereupon  it  is  easily  seen  that 


V  h  H  2rab  cos  <f>  * 


and  the  Doppler-shifted  frequency  is 


(3.6) 


(3.7) 


(3.8) 


£o[l  -  f  7'*Ab‘°S  <f>  -  f  RAB 


rnCOS 


However,  just  as  in  Case  I  we  may  drop  the  last  two  terms  in  the  bracket 
as  being  quite  negligible  in  comparison  with  the  first  term,  leaving 
simply 


f 


B 

n 


Proceeding  further  as  in  Case  I  we  find  the  signal  re-radiated  by 
scatterer  n  to  be ,  in  the  far  field  near  receiver  B, 
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->B  . 

E  -  o 
n 


Et  «P  tznlfa  C-  V  V'n  +  c<t'tn)])" 


with  t  being  the  nominal  time  of  arrival  of  the  re-radiated  wave  front 

from  0:  tat  +“(R.  +  U.  Here  a'  is  the  operator,  corresponding 
o  c  A  d 

to  (7  in  Case  X,  representing  symbolically  the  effect  of  the  scattering 
center  on  the  incident  pulse  for  re-radiation  in  the  primed  direction; 

CJ*  may  or  may  not  be  the  same  as  (7 ,  depending  on  the  nature  of  the 
scattering  center  involved.  Using  eq.  (3.8)  and  the  value  of  t  given 
above,  the  received  signal  B  becomes 


-B 

E 

n 


exp  [- 


4tt 


cos 


If  the  receiving  antenna  is  polarized  in  direction  [where  Rg  0] 

then  the  received  signal  from  scatterer  n  has  the  form 


n 


ei0o 

Vb 


Et  exp  [-  ^  iR 


AB 


cos  (-|)  ] 


(3.9) 


A  similar  analysis  can  of  course  be  done  for  remote  receiver  C  also; 
the  result  is  like  eq.  (3.9)  with  B  everywhere  replaced  by  C  and  with 
cos  P  =  R^*  Rc 

After  reception  the  return  signals  are  compared  with  a  reference 
signal  for  phase,  which  has  the  effect  of  eliminating  the  phase  factor 
0  but  inserting  a  "range-phase"  factor  which  has  the  form  4ttR^ /^q 
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in  the  monostatic  case  and  2tt(Ra+  Rg)/^  in  the  bistatic  case.  During 
the  course  of  the  data  processing,  however,  very  accurate  determinations 
of  range  are  made,  which  allow  the  range-phase  factor  to  be  removed  and 
also  permit  the  recorded  signal  amplitudes  to  be  compensated  for  the 
1/R|  or  l/RARg  dependence  thereby  transforming  the  results  into  absolute 
units.  In  the  remaining  sections  of  this  paper  it  will  be  taken  for 
granted  that  both  of  these  processing  steps  have  been  done. 
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SECTION  IV 


ISOTROPIC  POINT  SCATTERERS 

The  term  "isotropic  point  scatterer"  will  be  used  here  to  refer  to 
any  scattering  center,  whatever  its  nature,  which  satisfies  the  following 
cr  iter  ia . 

1.  The  scattering  center  can  be  represented  as  a  point  rigidly 

attached,  by  a  vector  r  (t)  ,  to  the  object’s  center  of  mass,  0,  and 

n 

rigidly  rotating  about  0  with  the  object's  angular  velocity  uj(t)  .  The 
magnitude,  r^ ,  of  the  vector  r^  is  a  constant,  independent  of  t. 

2.  In  the  linear  polarization,  the  scattered  wave  has  the  same 
polarization  as  the  incident  wave;  i.e.,  there  is  no  depolarization. 

This  is  assumed  to  be  true  both  for  the  monostatic  case  (back-scattering) 
and  for  the  bistatic  case  for  at  least  the  small  bistatic  angles  considered 
here . 

3.  The  scattering  is  characterized  by  a  complex  number,  the 

iY 

scattering  coefficient,  of  the  form  H^e  n  ,  where  gives  the  relative 

strength  of  the  re-radiated  signal  and  Y^  gives  the  phase  change  (if 

any)  upon  re-radiation.  The  numbers  H  and  Y  may  differ  from  one 

n  n 

scattering  center  to  another,  but  they  are  constant  in  time  and  angle. 

Suppose  now  that  a  unit  signal  is  transmitted  with  linear  polari- 

A  A 

zation  ET  =  e^  .  By  our  assumptions  above,  the  effect  of  scatterer  n 
on  this  signal  will  be  represented  (in  either  the  monostatic  or  the 
bistatic  case)  by 
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(4.1) 


iY 


CTE  =  o'E  -  e,H  e 
T  T  j  n 


Thus,  distinguishing  between  the  monostatic  and  the  bistatic  cases,  if 
the  receiving  antenna  has  linear  polarization  given  by 


(monostatic) : 


<  *  4i 


(bistatic) 


ER  ~  6i 


(4.2) 


then  the  received  signal  is  characterized  by 


(monostatic):  E^ 


iY 


R  *  CTET  "  ei  •  ej  Hne 


(4.3a) 


(bistatic) 


Ep  *  ct'E  =  e!  *  e.H  e  n 
k  l  1  j  n 


iY 


(4.3b) 


Making  use  now  of  Eq.  (3.5),  the  comment  at  the  end  of  section  III, 

the  definition  of  Scattering  Matrix,  and  the  fact  that  e.*  e.  =  6 . . ,  we 

1  J  ij 

see  that  in  the  linear  basis  and  for  the  monostatic  case  the  scattering 
matrix  for  the  nth  scatterer  has  the  elements 


iY  r  ^ 

(monostatic):  S  =  6.  .H  e  n  exp[-4iri  “  R.  *  r  ].  (4.5) 

ij  ij  n  A  n 

o 
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For  the  bistatic  case  the  situation  is  made  somewhat  more  complicated 


by  the  presence  of  the  matrix  e^  •  e^  which  is  not  simply  equal  to  6^ ^ . 
However,  for  small  bistatic  angles  it  can  be  shown  (see  Appendix  A)  that 


ej  “  6ij  "  (2)2  Kij 


(4.5) 


where  the  matrix  K  has  elements  of  order  unity.  Therefore  if  we  make 
use  of  Eq.  (3.9)  and  ignore  terms  of  order  (^)2  ,  we  see  that  in  the  linear 
basis  and  for  the  bistatic  case  the  scattering  matrix  for  the  n^  scatterer 
has  the  elements 


(bistatic):  S 


ij 


i  Y  r  ^ 

6^  Hne  exp[-4rri  ~  R 


?  cos 
n  z 


AB 


(4.6) 


Here  we  have  retained  the  factor  cos  (^)  in  the  exponent,  even  though 
it  is  approximately  equal  to  1  -  -|(-^ ) 2  ,  because  it  appears  multiplied 
by  what  may  be  a  large  number  (i.e.  4tt  r^/A^  and  phase  factors  are  always 
taken  modulo  2rr  so  this  may  make  a  significant  difference. 

Transforming  into  the  circular  basis  by  means  of  Eq.  (2.14),  we 
find  that  the  scattering  matrix  for  the  nth  scattering  center  has  the 
elements 


(monostatic)  : 


iY  r  „ 

cij“  (1  •  V  V  "  e*pt-4ni  RJ 


(4.7) 


iY  r 

(bistatic):  (1  -6ij)Hne  n  exp[-4ni  RAB*  rncos(^)].  (4.8) 
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These  values  of  the  scattering  matrix  elements  are  for  one  pulse 
only.  It  will  now  be  necessary  to  find  expressions  for  the  time-depend¬ 
ence  of  the  .  Since,  by  the  defining  properties  stated  previously, 

the  quantities  H^,  Yn>  and  r^  are  constants,  the  only  functions  which 
can  contain  the  time  dependence  of  C  are  R  •  in  Eq .  (4.7)  and 

~  O 

V  ?n  cos  (^)  in  Eq.  (4.8).  Let  us  examine  them  further  in  this 
regard . 

We  will  eventually  need  to  take  the  finite  Fourier  Transform  of 

the  elements  C_(t),  integrating  over  a  time  interval  of  length  T 

T  T 

centered  on  some  time  t:  t-^^t^t+-r  .  Usually  T  will  be 

o  o  2  o  2 

"small,"  the  precise  meaning  of  which  will  become  clear  as  we  proceed. 
Therefore  we  will  expand  all  time-dependent  functions  in  Taylor  series 

k 

about  t  and  keep  only  first-order  terms  in  small  quantities. 

°  dr^ 

For  r^(t)  we  note  that  — =  uu(t)  X  r^Ct)  >  whereupon  the  Taylor 
Series  gives 

rn(t)  ~r(t  )  +  (t-t  )  uu(t  )  X  r  (t  )  +  0  (u,  f)2  .  (4.9) 

n  no  o  o  no  2 

In  order  for  this  linearization  to  be  valid,  it  is  necessary  that 
T  a 

(uu  "2)  «  1  hold  true.  To  be  concrete,  suppose  we  require  that 

0)  T/2  <0.1  .  Then  writing  uu  =  2tt/t,  where  T  is  the  period  of  the 

object’s  rotation,  this  condition  becomes  approximately 

* 

For  certain  purposes  it  may  be  useful  to  retain  the  second-order  terms 
also,  but  this  will  not  be  considered  in  the  present  paper. 


30 


T  £  0.03  T  . 


(4.10) 


This  may  be  regarded  physically  as  a  limitation  on  the  length  of  the 
integration  interval  T  in  comparison  with  the  body  rotation  period  t, 
such  that  the  Fourier  Transform  not  be  "smeared-out"  by  rotation  through 
too  large  an  aspect  angle  during  the  integration.  If  amplitude  weighting 
is  used  during  integration  to  reduce  the  contributions  of  data  points 
near  the  ends  of  the  interval  (|t-to|~-|)  in  comparison  with  data  points 
near  the  center,  then  presumably  this  condition  may  be  relaxed  a  bit. 

How  much,  depends  on  the  form  of  weighting  used;  perhaps  T  <  0.05  t 
would  be  alright. 

A 

The  function  R^(t)  varies  quite  slowly  under  most  conditions;  if 
we  write 

dF  V  V*>  X  *A(t)' 

-2 

then  WA  <  10  rad/sec  except  near  zero-Doppler  of  very  close  passes. 
Furthermore,  UJA  itself  changes  quite  slowly,  so  that  it  is  reasonable 
to  assume  that  ^  «  0  for  tQ-  |  *  t  £  tQ  +  j  even  if  T  is  not 

absolutely  small.  Therefore  we  can  approximate  R  (t)  by  the  first 

n 

two  terms  of  its  Taylor  Series  expansion: 

RA(t)  -  RA(tQ)  +  (t-to)  wA(to)  X  RA(tQ)  +  0(wA  f)2  .  (4.11) 
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For  the  bistatic  case,  if  we  write 


dT  [*ABC0S  <2)]  S  “AB(t)  *  ®AB(t) 


and  make  use  of  Eq.  (3.8)  to  express  R.„  cos  (&)  in  terms  of  R.  and 

Ad  Z  A 

R,, »  then  it  turns  out  that  the  difference  between  tu,„  and  ul  is  of 
B  AB  A 

second  order  in  small  quantities.  Therefore  we  can  write,  analogously 
to  Eq.  (4.11)  , 


^  **AB(to^  +  X  *AB^  +  0('UJA  2^# 


(4.12) 


Substituting  Eqn.s  (4.9)  and  (4.11)  into  (4.7)  leads  to  the 
following  expressions  for  the  time  dependence  of  the  scattering  matrix 
elements  for  the  nth  scattering  center  in  the  monostatic  case: 


if  r 

Clj(t)  “(1-6ij>  V  ”  exp(-4TTl  J[iA(to)  •  fn(t0)-<t-t0)RA(t0)  x  fn(t0).n(t0>:j. 

(4.13) 


where  the  superscript  "M"  reminds  us  that  it  is  for  the  monostatic  case, 
the  exponential  has  been  linearized  in  terms  of  small  quantities,  and 


Q(t)  =  uu( t)  -  iyt) 


(4.14) 
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Similarly  for  the  bistatic  case  (with  superscript  "B"  as  indicator): 

B  if  r 

Cu  (t)  “  (1  -  6ij)H„e  n^Pt-4ni-j£  cos  (f)[RAB(to)-  rn(to) 

'  (t'to)“AB(to)  X  V'o*  '  S  (4-15> 

Suppose  now  that  the  target  object  is  made  up  of  N  such  isotropic 

point  scatterers,  and  that  as  ideal  points  they  never  block  or  shadow 

one  another.  Then  the  total  scattering  matrix  for  the  entire  object  is 

found  by  letting  each  element  be  the  sum  of  N  terms  of  the  form  (4.13) 

or  (4.15).  Or  alternatively  the  total  scattering  matrix  may  be  thought 

of  as  the  sun  of  N  scattering  matrices,  each  with  elements  of  the  form 

(4.13)  or  (4.15).  These  two  pictures  are  mathematically  equivalent, 

and  represent  in  terms  of  this  simple  scattering  center  model  what  one 

should  expect  to  measure  as  the  scattering  matrix  over  the  time  interval 

from  t  -  to  t  +  ^  . 

o  Z  o  2 

For  the  measured  scattering  matrix  we  attempt  to  resolve  the 

individual  scattering  centers  in  "Doppler"  by  performing  a  finite  Fourier 

Transform  of  the  matrix  elements.  Let's  see  what  this  corresponds  to 

in  terms  of  the  model.  The  finite  Fourier  Transform  of  a  function  S(t) 

on  the  interval  t  -  t  i  t  +  |  is  defined  to  be 
o  z  o  Z 
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S(f) 


S(t)  e-2TTif(t‘to) 


dt 


(4.16) 


Operating  thus  on  the  sum  of  terms  of  the  form  (4.13)  we  find  that 


N  .  A 

10, 


^j( f)  =“  (l-6i:j)  T  ^  Hne  n  [TT  T(f'fn^  *  <4-17> 


n=l 


where 


H  "  4tt  7^  Ra(C  )  *  r  (t  )  , 
n  n  Ao  A  o  nv  o' 


(4.18) 


a  2r  - 

C  =  l"2  R.(t  )  X  i  (t  )  •  3(t  ) , 

n  ao  A  o  n  o  o 


(4.19) 


and  j  00  t^ie  Spherical  Bessel  Function  of  order  zero,  given  by 


jQ00  =  -X  s  sine  (x/n)  .  (4.20) 
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The  equivalent  expression  in  the  bistatic  case  is 


N 

i(A 


C^(f)  -  <1-6±J)  Hne  n  jQ[n  T(f-f®)]  ,  (4.21) 


n*l 


where 


0B  =  Y  -  4n^  R.R(t  )  *  r  (t  )  cos  (|) 
n  n  A0ABo  no  / 


(4.22) 


and 


fB  h  2  y1  R  (t  )  X  r  (t  )  *  0  (t  )  cos  (|). 
n  A  AB  o  n  o  o  ^ 

o 


(4.23) 


The  Fourier  Transformed  scattering  matrix  elements  (f)  corres¬ 

pond  to  the  results  of  the  "Doppler  mapping"  process.  They  will  be  dis¬ 
cussed  in  some  detail  in  a  later  section,  but  at  the  present  time,  we 
wish  to  comment  on  just  one  thing:  the  fact  that  each  term  of  the  sums 
in  (4.17)  and  (4.21)  contains  a  function  of  the  form  F^(f)=  THnjo[nT(f-fn) ] . 
This  function  has  a  strong  main  peak  at  f-fn  and  lesser  subsidiary  peaks 
in  the  sidelobes;  the  amplitude  of  the  main  peak  is  given  by  TH^  and 
its  width  is  proportional  to  (TT  T)  .  This  behavior  is  the  explanation 
for  the  fact  that  the  modulus  of  the  "Doppler  maps"  contains  sharp  peaks 
at  values  of  the  frequency  corresponding  to  individual  resolved  scattering 
centers . 
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A  bit  of  further  realism  can  be  added  to  this  model  if  we  allow 
for  the  possibility  of  shadowing.  In  a  real  object,  made  up  of  non¬ 
ideal  scatterers,  there  is  a  good  chance  that  during  the  time  of  integra- 

tion  1 1 - 1  I  £  -r  some  of  the  scattering  centers  will  not  be  seen  at  all 
1  o  2. 

(being  "in  shadow'1  the  whole  time)  and  others  will  be  seen  only  part  of 

the  time  (entering  or  emerging  from  the  shadow  at  some  points  in  the 

integration  interval).  The  completely  shadowed  point  is  easy  to  treat: 

it  simply  makes  no  contribution  at  all  to  the  scattering  matrix.  What 

about  the  other  case?  Suppose  a  point  "V%  is  visible  from  t-t  -  until 

T  T 

t=t  +  t0,  where  --r^-t.  <t«^T,  and  is  invisible  during  the  rest 
O  Z  Z  I  z  z 

of  the  interval.  In  this  case  a  straightforward  integration  shows  .that 
the  transformed  scattering  matrix  has  the  form 


*0# 


C_(f)=“  (l-6„)  expC-inCtg-tj)  (f-f^)  ]jo[nX£  T(f-f^)],  (4.24) 


where  =  (tx+  t3)/T  represents  the  fraction  of  the  total  integration 

A  B 

interval  for  which  the  scatterer  is  visible,  and  where  0^  =  0^  or  0^ 
and  f^  =  f^  or  f®  depending  on  whether  the  situation  is  monostatic  or 
bistatic . 

Clearly,  then,  the  effect  of  the  partial  shadowing  is  (a)  to  decrease 
the  amplitude  and  increase  the  width  of  the  corresponding  peak  of  the 
Doppler  map  modulus  in  proportion  to  the  fraction  of  the  integration  time 
spent  in  shadow,  and  (b)  to  insert  a  phase  change  proportional  to  ( t t vf-f^) > 
so  that  there  is  no  phase  change  if  the  visible  "window"  is  evenly  centered 
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in  the  integration  interval,  and  in  any  case  there  is  no  phase  change 
at  f=f^:  the  center  of  the  peak. 

Further  comments  concerning  these  scattering  matrix  elements  will 
be  reserved  for  a  later  section. 
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SECTION  V 


FIXED  POINT  DIPOLE  SCATTERERS 

The  simplest  model  of  a  scattering  center  which  depolarizes  and 
has  "orientability"  is  the  fixed  point  dipole,  by  which  we  mean  a  point 
dipole  of  dipole  moment  p^(t)  rigidly  connected,  by  a  vector  r^(t)  ,  to 
the  object’s  center  of  mass,  0,  and  rigidly  rotating  about  0  with  the 
object’s  angular  velocity  uu(t)  .  The  magnitudes  r^  and  p^  are  constant, 
independent  of  t.  When  we  say  that  the  vector  p^  is  rigidly  connected 

“♦  —4  “♦ 

by  r  to  0,  we  mean  that  p  transforms  in  the  same  way  as  r  under 

rn  y  n 

rotation  about  0,  i.e. 


dt  Pn  =  X  (5.1) 

This  may  either  be  regarded  as  a  definition  of  ’’fixed  point  dipole",  or 
may  be  shown  to  follow  necessarily  from  a  mathematical  formulation  of 
what  one  might  intuitively  mean  by  the  statement  that  p^  is  "rigidly" 
connected  by  r^  to  0:  See  Appendix  B. 

Suppose  now  that  a  unit  signal  with  linear  polarization  e^ 

is  transmitted  by  antenna  A  toward  the  target.  The  excitation  of  dipole 
n  by  this  signal  will  be  proportional  to  E^*  p^.  The  dipole  then  begins 
to  radiate,  aid  the  appropriate  formula  for  the  far-field  radiation  of  an 
oscillating  electric  dipole  tells  us  that  as  it  approaches  the  receivers 
the  reradiated  signal  is  proportional  to  the  excitation  and  to 
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in  the  monostatic  case 


RA  X  (pn  x  V 


and  to 


A  — ♦ 

Rfi  x  (pn  x  in  the  bistatic  case, 


Therefore,  in  the  notation  of  Section  III, 


0ET  '  (eY  fn>[*A  x  (pn  x  VI 


(5.2a) 


and 


°,EIS  (Y  Y[V  <?„  X  V] 


(5.2b) 


If  the  receiving  antenna  has  linear  polarization  given  by 


(monostatic)  : 


(bistatic)  :  RR  =  *  (5.3) 

and  we  recall  that  e^'  RA  =  0  and  e^  *  Rg  =  0,  then  the  received 
signal  is  characterized  by 

(monostatic):  •  CTET  *  pn^®j"  Pn^®i’  Pn^  (5.4a) 

(bistatic)  :  E®  '  ct'Et  =  p^e^*  Pn)(e^  *  Pn)  •  (5.4b) 
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Making  use  now  of  eq.  (3.5),  the  comment  at  the  end  of  Section  III, 
and  the  definition  of  Scattering  Matrix,  we  see  that  in  the  linear 
basis  and  for  the  monostatic  case  the  scattering  matrix  for  the  nth 
dipole  has  the  elements 


Pn(V  Pn)(V  pn>  el,P[‘4nl  Y  *„]• 

J  n 


(5.5) 


For  the  bistatic  case  the  scattering  matrix  elements  are 


Pn(Si'  Pn)(V  p„>  f  Pn  cos  <!> 


(5.6) 


Notice  that,  as  required, 

*aiT*a  as  0 

Now  we  must  make  the 
define 


Lim  SB.  =  SM. 
@-0  1J  1J 


transformation  into 


since  e!  -  e.  and 
i  i 


the  circular  basis. 


We 


(5.7) 


and  use  eqns  (2.14)  to  find  the  scattering  matrix  elements  in  the  monostatic 
case : 
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(5.8a) 


,M 


-it 


M 

C11  =  *Pne  n  C(V  pn)2  "  (V  Pn)2+  2i(V  Pn)(V  Pn> 3 


— M  _M  i  "n  r  / A  A  v  Z  .  a  a  .  2  n 

C12  =  C21  ■  ^„e  [<V  Pn>  +  <V  p„)  ] 


-*e 


-  v2 


(5.8b) 


C22  =  ^n6  n  C(V  Pn)2_  (V  Pn)2_  2i(V  Pn>  ( V  Pn)]*  (5‘8c) 


The  bistatic  case  is  somewhat  more  complicated:  because  of  the 
presence  of  the  primes  in  eq.  (5.6)  the  expressions  found  upon  using  eqns . (2 . 14) 

T>  -n 

do  not  simplify  very  much  since,  for  example,  S12  T  S21  .  But  for  small 

bistatic  angles  these  difficulties  can  be  removed.  Suppose  we  define  basis 

vectors  e^’  in  a  plane  normal  to  R^g  in  analogy  to  the  definitions  of 

e^  in  a  plane  normal  to  R  .  Then  it  can  be  shown  [see  Appendix  A]  that 
B  2 

to  order  0(p  the  elements  of  the  scattering  matrix  in  the  circular  basis 
and  the  bistatic  case  can  be  found  from  those  in  the  monostatic  case  by 
simply  replacing  e^  by  e£  and  by  where 

?n  *  4"  f  “s  <!>  •  <5‘9> 

o 

This  statement  may  be  considered  to  be  the  form  of  the  monostatic-bistatic 
equivalence  theorem  for  a  point  dipole  scatterer.  It  should  be  carefully 
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noted  that  this  equivalence  does  not  hold  true  for  the  scattering  matrix 
in  the  linear  basis. 

Hereafter  we  shall  consider  only  the  raonostatic  case,  the  bistatic 
one  being  deriveable  from  it  by  the  above  transformations. 

The  scattering  matrix  elements  can  be  expressed  in  a  simple  way  in 

* 

terms  of  meaningful  geometrical  quantities.  The  plane  defined  by  e^ 
and  e2  is  often  called  the  radar  polarization  plane.  The  projection  of 
the  vector  p^  onto  this  plane  makes  an  angle  9n  with  the  e^-axis,  and 
this  angle  is  often  called  the  Polarization  Angle  of  the  dipole.  Similarly, 
the  projection  of  p^  onto  the  vector  RA  (normal  to  the  Polarization  Plane) 
defines  the  Aspect  Angle,  £n>  of  the  dipole.  See  Figure  3  below: 

A 

R* 


Figure  3.  Decomposition  of  the  Vector  Pn  . 
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Notice  that 


cos  Q 

n 


e.  *  p  =  sin  C  cos  0 
1  rn  n  n 


V  '  sin  Cn  sin  9, 


(5.10) 


If  we  substitute  eqns  (5.10)  into  (5.8)  we  easily  find  the  scattering  matrix 
elements 

-(§  -  20  ) 

CH  -  3i"  Cne 

C12  *  C21  =  Bia\  e 

,  -i(?  +  26  ) 

C22  =  Sin  ^n6  *  (5.11) 

These  can  all  be  stated  in  one  equation: 

Cij  =  ^Pn  Sin2  exPt"1[?n+  26_(-l)j  6n  ^  •  (5.12) 
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The  next  step  is  to  find  the  time  dependence  of  these  matrix  elements. 
Now,  the  time  dependence  of  has  already  been  discussed  in  Section  IV 

[see  eqns  (4.13)  and  (4.15)].  The  factor  sin  <J^(t)  can  be  expanded  in 
Taylor  series,  making  use  of  eqns  (5.10),  (4.10),  and  (5.1)  along  the 
way,  and  linearized  in  terms  of  small  quantities  as  has  been  done  before; 
the  result  is 


Sin2<»n(t)  “■  sin2^n(to)  +  2(t-to)  COS^n(to) 


MO 

A  o 


Pn(t  ) 

n  o 


**„>• 


(5.13) 


and  finally  we  may  write  for  8n ( t )  the  expression 


e  (t)  -  e  (t  )  +  (t-t  )  e  (t  ) 

n  n  o  o  n  o 


(5.14) 


That  the  linearization  here  is  legitimate  may  be  shown  by  expressing 
9^  in  terms  of  more  fundamental  quantities;  it  turns  out  to  be  of  first 
order  in  small  quantities  like  u)  and  uu^,  with  higher  derivitives 
being  of  higher  order.  The  complete  expression  for  9^  is  not  simple, 
but  it  may  be  important  in  later  stages  of  analysis  to  know  its  form, 
so  it  is  derived  and  discussed  in  Appendix  C. 


Now  the  linearized  time  dependence  of  C_  can  be  written  down. 
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Let 


f  -  2  r1  Mt  )  x  r  (t  )  *  fi(t  ), 
n  A  A  o  no  o 

o 


(5.15) 


gnn(to)  sin  Cn(to)  =  RA(tQ)  X  pn(to)  *  fi(to), 


(5.16) 


and 


K  * 4n  r  ftACto)  '  cos  <f> 

o 


(5.17) 


A  A 


where  R.  *-  R.  _  in  the  bistatic  case 
A  AB 


and  @  =  0  in  the  monostatic  case. 


Then 

Cij(t)  “*  ^Pnsin^(to)Csin^n(to)  +  2(t-to)8nn(to)  cos^n(to)] 
exp{-i[0n+  26i.(-l)jen(to)]}  exp{2TTi(t-to)[fn-  ^  6^ (-l)j  6n(tQ)  ]}  .  (5.18) 


Suppose  now  that  the  target  is  made  up  of  M  such  point  dipole  scatterers, 
and  that  as  ideal  points  they  never  block  or  shadow  one  another.  Then  the 
total  scattering  matrix  for  the  entire  object  is  found  by  letting  each 
element  be  the  sum  of  M  terms  of  the  form  (5.18).  In  terms  of  our  simple 
scattering  center  model  this  represents  what  one  should  expect  to  measure 
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T 

as  the  scattering  matrix  over  the  time  interval  |t-to|s  — 

Applying  the  finite  Fourier  transform  (4.16)  to  this  sum  of  terms 
we  find  that 

C .  (f)  Gm(f)  6  (5.19) 

m=l 

where 


G  (f)  =  ^Tp  sin  f  (t  )  {sin£(t  )  j  [nT(f-f*)  3 
m  m  Tn  o  to  o  o  m 

-  iTgmn(to)  cos^(to)j1[Trr(f-f *)]),  (5.20) 

with 

0*  =  0  +  26.  .(-i)j  e  (t  ), 

m  *m  ij '  '  mv  o’ 

f*  =  f  .  I  6  .(-1)J  0  (t  ), 
m  m  TT  ij  mv  o7> 

and  j^OO  the  Spherical  Bessel  Function  of  order  one,  given  by 

J]_GO  =x  ~  cos  x)*  (5.23) 

These  Fourier-Transformed  scattering  matrix  elements  C^j(f)  should 
correspond  to  the  results  of  the  Doppler  mapping  process.  They  will  be 
discussed  in  some  detail  in  a  later  section,  but  at  the  present  time  we 


(5.21) 

(5.22) 
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wish  merely  to  point  out  that  the  modulus  of  Gm(f)  usually  has  a  peak  at 
★ 

f*f  ,  though  both  the  magnitude  and  the  width  of  this  peak  vary  strongly 

with  the  aspect  angle  Q  (t  )  and  for  sufficiently  small  values  of  C 

mo  m 

it  may  be  so  small  and  broad  as  to  escape  notice.  The  rather  complicated 

function  G  (f)  will  be  examined  more  fully  later, 
m 

•k 

Notice  that  when  a  recognizable  peak  does  occur  at  f=f  for  a 

m 

particular  scatterer,  it  does  not  appear  in  the  same  place  in  each  of 
the  four  Doppler  maps.  its  location  will  be  the  same  for  the  two  off- 
diagonal  elements,  and  it  will  appear  in  the  diagonal  elements  shifted 
by  equal  amounts  in  opposite  directions.  This  shift  may  make  the  process 
of  associating  corresponding  peaks  more  difficult,  but  on  the  other  hand 
if  the  association  can  be  made  then  the  shift  provides  useful  information 
about  the  dipole  involved: 


Af  =  f  -  fn  =  ;  6  (t  ). 
m  22  11  TT  m  o 

Similarly  the  shift  in  the  phases  corresponding  to  these  two  diagonal 
peaks,  if  it  can  be  determined,  yields  the  polarization  angle  directly: 


Hi  5  *22*  0U=  •40n(to)  • 

As  in  the  case  of  isotropic  point  scatterers,  further  realism  can  be 

added  to  the  model  by  allowing  for  the  possibility  of  shadowing.  Once 

again,  as  at  the  end  of  Section  IV,  we  consider  a  scattering  center  "V1 

T 

which  is  visible  from  t=t  -  t-  till  t=t  +  t0,  where  -  “  ^  -t-<  t0  ^ 

o  I  o  z  z  i  z 
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HI  CM 


and  t^+  =  X^T-  Then  the  transformed  scattering  matrix  for  this 

scatterer  takes  the  form 

★ 

~  -10.  * 

Ctj(f)  =“  %y^Tp^sin  C^(tQ)e  exp[-iTT(t2-  {[sin  C^(tQ) 

+  (t2-  tx)  &l  n(to)  cos  c4(te)]  i0^T(f-f*)] 

-  ix^Tg^  n(to)  cos  Qt(t0)  j1[TTx>tT(f-f*)]}  .  (5.24) 


The  same  comments  which  were  made  (at  the  end  of  Section  IV)  about  the 
results  of  shadowing  for  isotropic  point  scatterers,  Eq.  (4.24),  also 
hold  true  about  the  above  results  of  shadowing  for  point  dipole  scatterers. 
The  modulus  peaks  are  reduced  in  magnitude  and  increased  in  width,  and 
a  phase  term  is  inserted  which  vanishes  when  f=f^  and  which  is  pro¬ 
portional  to  the  "asymmetry  of  the  shadowing":  t2~  t^. 
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SECTION  VI 


SOME  REMARKS  CONCERNING  THE  DOPPLER  MAP  FUNCTIONS 


Consider  now  a  complex  target  composed  of  N  scatterers,  some 
being  isotropic  points  as  described  in  Section  IV  and  the  others  being 
point  dipoles  as  described  in  Section  V.  The  transformed  scattering  matrix 
or  "Doppler  Map  Function"  for  this  target  may  be  written  in  the  form 


exp  j  iAn(ij,  f)  | 


(6.1) 


where  the  functions  M  and  A  depend  in  form  on  whether  scatterer  n 

n  n 

is  an  isotropic  point  or  a  point  dipole,  and  where  the  bistatic  cases  are 
referred  back  to  the  appropriate  monostatic  situations  by  invoking  the 
monostatic-bistatic  equivalence  theorem.  Allowing  for  the  effects  of 
shadowing,  these  functions  are: 

(1)  for  isotropic  point  scatterers,  using  eq.  (4.24): 


(6.2) 


An(ij,f)  =  0n-  n(t2-t1)(f-fn)  , 


(6.3) 


with 


2r 

n 


(6.4) 
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(6.5) 


0n  ■  V  4"  f  VCo>  '  fn(eo)  cos<§>  ■ 

o 


where  *■  ft^g  in  the  bistatic  case  and  @  =  0  in  the  monostatic  case; 


(2)  for  point  dipole  scatterers,  using  eq.  (5.24): 


Mn(ij,f)  “  %  XnTpn  sin  £n( tQ)  |Gn(ij,f)|  (6.6) 

An(ij,f)  =  0*  (ij)  -  TT(t2-t1)[f-f*(ij)  ]  +  Yn(ij,f),  (6.7) 


|Gn(ij,f)|2  =  [sin  Cn(tQ)  +  (t2-tl)  gnO(to)  cos  Cn(to)]{jo[nXnT(f-f*(ij))]}: 


+  CXnTgn  n(to)  cos  Gn(to)]2{j1[nXnT(f-f*(ij))]}2.  (6.8) 


tan  Yn(ij,f) 


XnTgn  n<0«>«  G_(t  ) 
n  n _ o _ n  o 


j1CTTXnT(f-fn(ij))] 


(6.9) 


fn(iJ>  *  f„  -  9  Vl)J  MO, 


n  o 


(6.10) 
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; 


o 


where 


AB 


in  the  bistatic  case, 


and 


M  0  in  the  monostatic  case 


Equation  (6,12)  defines  the  quantity  gn,  whose  absolute  value  is  less 
than  unity. 

It  has  been  found  most  convenient  to  display  experimentally- 
determined  "Doppler  Maps**  in  the  form  of  modulus  and  argument  of  the 
complex  function  plotted  against  frequency  f.  So  for  comparison  with 
experiment  we  express  eq.  (6.1)  in  the  form 


C^Cf)  =-Mij(f)  exp  {LA^Cf)} 


(6.13) 


where  it  is  not  difficult  to  see  that 


N  N 


(6.14) 


n~l  m“l 


and 


tan  A±j (f) 


Mn(ij,f)sin  An(ij,f) 
N 

Mn(ij,f)cos  An(ij,f) 

n»l 


(6.15) 


The  behavior  of  M^(f)  and  A^(f)  as  functions  of  f  depend 

primarily  on  the  form  of  the  functions  M^Cijjf).  Qualitatively,  as 

has  already  been  remarked,  these  functions  usually  (though  not  always: 

see  below)  are  more  or  less  sharply  peaked  at  f*=f^  and  relatively 

small  elsewhere;  i.e.,  they  are  rough  approximations  to  6(f-f  ).  Even 

n 

if  for  certain  types  of  scattering  centers  the  functions  Mn(ij,f) 
do  not  have  the  precise  forms  shown  in  Equations  (6.2)  or  (6.6),  it  is 
reasonable  to  assume  that  in  most  cases  they  will  have  the  same 
qualitative  properties,  so  that  the  argument  below  is  still  applicable. 

Now  if  Mn(ij,f)  were  indeed  very  close  to  a  true  delta  function 
(which  would  be  the  ideal  situation) ,  what  should  we  expect  to  happen? 
Well,  equations  (6.14)  and  (6.15)  would  then  be  effectively  "decoupled" , 
in  the  sense  that  we  would  have 


Vf) 


M  (ij,f  )  at  f“f  for  any  n 
n'  J’  n'  n  } 

~  0  at  f  /  f  for  any  n 
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A  (ij,f  )  at  f=f  for  any  n 
n  n  n  J 

random  at  f^f^  for  any  n 

Thereupon  by  observing  the  locations  of  the  peaks  in  the  modulus  of  the 
Doppler  Map  function  we  could  easily  find  the  true  values  of  f^  and  0^ 

(mod  2n)  and  be  well  on  the  way  towards  solving  for  the  locations,  motions, 
and  orientations  (if  applicable)  of  all  the  scatterers.  [if  scatterer 

“fc 

a  is  a  dipole,  substitute  f  and  0  for  f  and  0  everywhere  in 
this  paragraph.] 

But  if  Hn(ij,f)  is  not  very  close  to  a  perfect  delta  function, 
but  rather  has  significant  sidelobes,  then  the  different  terms  of 
equations  (6.14)  and  (6.15)  become  seriously  cross-coupled.  This  inter¬ 
modulation,  as  it  may  be  called,  has  three  deleterious  effects,  all  of 
which  have  been  easily  observed  in  simulated  data  made  up  to  satisfy 
Eq.  (6.2).  These  are: 

1.  Spurious  peaks  may  appear  in  the  modulus  function  at  values 
of  f  not  corresponding  to  any  true  scatterer. 

2.  The  values  of  f  at  which  genuine  peaks  appear  in  the  modulus 

★ 

function  are  somewhat  displaced  from  the  correct  values  f  or  fR 

3.  The  values  of  the  argument  corresponding  to  peaks  in  the 

•k 

modulus  function  are  somewhat  displaced  from  the  correct  values  0  or  0  . 

n  n 
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These  problems  can,  if  sufficiently  severe,  make  correct  inter¬ 
pretation  of  the  data  rather  difficult.  Therefore  it  is  important  to  try 

to  ensure  that  the  functions  M  (ij,f)  are  reasonably  good  approximations 

n 

to  delta  functions.  Naturally,  we  can't  change  the  data  as  the  target 
presents  them  to  us.  But  it  may  be  possible  to  take  steps  during  the 
Fourier  integration  procedure  to  enable  the  functions  to  have  peaks 

as  narrow  and  sidelobes  as  low  as  the  natures  of  the  individual  scattering 
centers  allow.  Amplitude  weighting,  for  example,  can  considerably 
reduce  the  sidelobes,  although  at  the  cost  of  some  broadening  of  the  peaks. 
Increasing  the  integration  time  T  with  the  help  of  previously  estimated 
body  motion  constraints  might  go  a  long  way  toward  alleviating  this 
problem,  if  it  can  be  done. 

After  the  above  discussipn  involving  the  general  porperties  of  the 

functions  Mn(ij,f),  this  seems  a  good  place  to  illustrate  the  specific 

forms  of  the  moduli  | |  for  certain  typical  cases. 

First  of  all,  for  isotropic  point  scatterers  Eq.  (6.2)  shows 

that  |Mn(ij,f)|e  0  for  ij=*ll  or  22,  and  that  |Mn(ij,f)|*  TXnHnjJ}nTx  n(f-fR)  ] 

for  ij  *  12  or  21.  Figure  4  shows  a  graph  of  |jQ(x)|  and  |j^(x)|:  it 

is  easy  to  see  that  the  maximum  value  of  | |  is  given  by  TxnHn  and 

4  50 

occurs  at  f“f  .  The  first  sidelobe  maximum  occurs  at  f  ^  f  +  —  , 

n  rr-  ttTv 

n 

and  its  value  is  approximately  given  by  0.22  TXnHn>  down  only  6.6  db 
from  the  peak.  The  3  db  width  of  the  main  peak  is  Af  ■  3.8/nTXn.  Of 
course  appropriate  weighting  reduces  the  sidelobe  level  considerably 
whilst  increasing  the  width  of  the  peak,  so  there  isn't  really  much 
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(  ABSOLUTE  UNITS  ) 


CD 


X  (  RADIANS  ) 

Figure  4.  ABSOLUTE  VALUES  OF  THE  FIRST  TWO  SPHERICAL 

BESSEL  FUNCTIONS 
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point  in  using  the  exact  values  just  quoted  above  except  as  indications 

of  how  the  moduli  | |  behave  as  functions  of  the  parameters  T,Xn  and 

H  . 
n 

For  point  dipole  scatterers  the  situation  is  considerably  more 
complicated  and  is  made  even  more  so  when  weighting  is  used.  To  get 
at  least  some  feel  for  the  behavior  of  |Mn(ij,f)|  under  several  different 
sets  of  parameters,  consider  Figures  5  through  9.  The  functions  plotted 
here  are 


|m(x)  |  =  I  sin  Q  |  [sin^Cj^Cx)  +  e2cos\ j^(x)  (6.16) 

for  various  values  of  the  parameters  Q  and  e,  using  the  same  relative 
scale  of  units  for  all  five  graphs.  Referring  to  Eqns.(6.6)  and  (6.8), 
it  is  clear  that  these  functions  represent  | [  for  fixed  values  of 
p^,  T,  and  y^,  plotted  against  the  independent  variable 

x » >v<f-0  • 

with  t^=  t^  and  e  =  xnTgnH.  Under  ordinary  conditions  the  parameter 
e  will  have  absolute  value  less  than  unity,  so  this  is  the  range  of  values 
used  in  Figures  5-9;  the  values  of  aspect  angle  Q  used  in  the  figures 
are  a  representative  sample:  10°,  30°,  45°,  60°,  and  80°.  Based  on 
these  figures  several  observations  should  be  made. 

(1)  Except  for  the  smallest  aspect  angles,  the  functions  all 
behave  qualitatively  like  | j  (x) |  in  that  they  have  a  main  peak  at  x*0 
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Figure  5.  THE  FUNCTION  lM  (X)  I  FOR  £  =  10°  AND 
VARIOUS  VALUES  OF  « 
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6.  THE  FUNCTION  lM(X)|FOR  £  «  30°  AND 
VARIOUS  VALUES  OF  c 
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50 


Figure  7.  THE  FUNCTION  I  M  ( X)  |  FOR  £=  45°  AND  VARIOUS  VALUES  OF  « 
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I.  THE  FUNCTION  IM(X)IfOR  C  ■  60#  AND 
VARIOUS  VALUES  OF  * 
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Figure  9.  THE  FUNCTION  |  M  (X )  |  FOR  £=  80*  AND  VARIOUS  VALUES  OF  < 
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and  smaller  sidelobe  peaks  at  or  beyond  x“4.5.  For  small  aspect  angles 
(e.g.  Q  =  10°)  the  main  peak  may,  if  e  is  large  enough,  occur  at  a  value 
of  x  greater  than  0. 

(2)  The  maximum  value  of  the  main  peak  increases  rapidly  with 
increasing  aspect  angle  Q.  Therefore  in  practice  the  anomalistic  peaks 
which  for  small  Q  occur  at  x>0  will  not  be  very  observable  in  comparison 
with  the  much  larger  and  more  "normal"  peaks  representing  dipoles  with 
larger  values  of  C,. 

(3)  For  a  given  value  of  aspect  angle  £  the  width  of  the  main 
peak  increases  with  increasing  e,  but  the  difference  becomes  smaller 
as  Q  becomes  larger,  so  that  when  C,  m  80°  the  function  |m|  is  almost 
independent  of  e. 
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appendix  a 


THE  MONOSTATIC-BISTATIC  EQUIVALENCE  THEOREMS 

Referring  back  to  Figure  1,  we  begin  by  defining  two  new  unit 
vectors,  S  and  t  : 


k  +  k' 

|k  +  k' | 

A  A 

k  -  k' 

|k  -  k*| 

from  which  it  is  not  hard  to  see  that 

k  +  k*  =  2S  cos  (—) 

A  A  /V  Q 

and  k  -  k'  =  2t  sin  (-*) 


S  = 


t  = 


(A-l) 


(A-2) 


Notice  that  S  is  a  unit  vector  in  the  direction  of  the  bisector  of 
the  bistatic  angle,  equivalent  to  R^g  of  Eq.(3.6)  and  others,  and 

A  A  A 

t  is  a  unit  vector  normal  to  S.  The  line-of-sight  vectors  k  and 
k'  can  be  expressed  in  terms  of  S  and  t  : 


k  i  S  cos  (y)  +  t  sin  (|-) 
k'=  S  cos  C^)  -  t  sin  Cj) 


(A-3) 
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1 


We  shall  also  assume  that  the  bistatic  angle  is  M small,"  meaning  that 

A  A 

we  can  neglect  all  but  the  lowest  powers  of  3/2  ,  and  that  n  *  n1 
to  sufficient  accuracy. 

Now  after  a  great  deal  of  inelegant  and  messy  vector  manipulations, 
starting  with  the  definitions  of  e^  and  ej  ,  it  turns  out  to  be  possible 
to  write 

e  -e"  +  uf+  0(f) 

1  1  2  2 

e  *  e  -  0  ,  +  0(f) 

1  1  2  2 

(A-4) 


e  =  e"  +  V  f  +  0(f)2 

2  2  2  2 

e'  =  e"  -  V  f  +  0(f)2  , 

2  2  ^  ^ 


where  the  vectors  e" 

i 

to  the  definitions  of 
where 


are  defined  relative  to  S  in  complete  analogy 

/V  A 

relative  to  k  ,  and 


U 


n»t  (n»S  )  (n*t)  e'j 
|n*s|  |n*S|2 


V  =  txe 


..  _  lai§l  t  + 


n*S 


(n*S)  (n*t)  e'' 
|nxs|2 


(A-5) 
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From  this  result,  it  is  obvious  that 


(A-6) 


K  = 


-*  ->  \ 
U-V  ' 


V 


(A-7) 


this  result  is  given  in  the  text  as  Eq.(4.5),  where  it  constitutes  the 
form  of  the  monostatic-bistatic  equivalence  theorem  for  an  isotropic 
point  scatterer. 

For  a  point  dipole  scatterer  the  situation  is  slightly  more  compli- 


cannot  be  done:  terms  of  order  6/2  remain  in  the  expressions  and  make 
the  equivalence  impossible.  This  means  that,  in  general,  the  monostatic- 
bistatic  equivalence  theorem  does  not  hold  true  for  orientable  scatterers 
in  the  linear  basis.  But  if  we  then  use  Eqns.(2.14)  to  transform  to  the 
circular  basis,  these  inconvenient  linear  terms  drop  out  of  all  the 

expressions,  leaving  the  equivalent  form  with  e”  replacing  e^  and 

B  M  B  ^ 

replacing  ,  good  up  to  a  term  of  order  (“)  .  Thus,  the  monostatic- 

bistatic  equivalence  theorem  does  indeed  hold  true  for  orientable  scatterers 

in  the  circular  basis. 
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APPENDIX  B 


DEFINITION  OF  RIGIDLY  CONNECTED  POINT  DIPOLE 


Consider  vector  p  rigidly  connected  to  the  end  of  vector  r  , 
both  vectors  having  constant  magnitude,  and  vector  r  rotating  about 
the  origin  with  angular  velocity  uj  .  We  know  that 


dr 

dt 


=  U) 


*  r  . 


(B-l) 


When  we  say  that  vector  p  is  rigidly  connected  to  vector  r  , 

— ♦ 

we  intuitively  mean  that  the  relationship  between  r  and  p  must 
remain  invariant  as  the  vectors  change  direction  with  time.  This 
means,  in  the  first  place,  that  the  angle  between  the  two  vectors  must 
not  change,  i.e. 


r(t)  *  p(t)  "  constant 


(B-2) 


And  in  the  second  place  it  means  that  the  vector  p  must  not  "twist" 

— ♦  — ♦ 

around  vector  r  as  they  change  with  time:  the  plane  defined  by  r 

and  p  should  not  rotate  around  r  .  This  condition  is  a  little  more 

difficult  to  express  mathematically,  but  we  shall  do  so  in  a  moment. 

The  condition  p  •  p  *  constant  leads  directly  to  the  requirement 
— ♦  — ♦ 

that  P  *  =  0  >  which  means  that  ^  can  always  be  written  in  the 

form 
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(B-3) 


dp  _  -  - 

-rf  =  U)  *  p 
dt  p  r 


for  some  vector  uu  .  Now  the  "no  twist"  condition  stated  above  can  be 
P 

put  mathematically  by  saying  that  uu  shall  have  no  component  along 


r  ,  i.e. 


cu  •  r  =  0 
P 


(B-4) 


By  differentiating  Eq.  (B-2)  and  making  use  of  Eq.  (B-l)  we  find 

that 


r*(^-<ju*p)=°  , 


which  implies  that 


= 


dt 


—  U)  x  p  +  a  x  r 


(B-5) 


for  some  vector  a  as  yet  unknown.  Taking  the  dot  product  with  p 

—4 

and  recalling  that  p  •  ^  =  0  we  see  that 


a  •  r  x  p  =  0  , 


which  says  that  a  lies  in  the  plane  of  r  and  p  ,  i.e. 


a  =  ar  -  bp 


for  some  numbers  a  and  b  .  Substituting  this  form  of  a  into 


67 


Eq.  (B-5)  we  get 


i 2 

dt 


=  (  uu  +  br  ) 


x  p 


which,  with  definition  (B-3),  indicates  that 


uj  =  id  +  br  . 

P 


When  this  form  is  tested  according  to  condition  (B-4)  it  becomes 
necessary  to  set  b  =  0  ,  whence 


ID 


P 


or 


dj 

dt 


OU  X  p 


(B-6) 


This  is  the  desired  result. 


« 
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APPENDIX  C 


THE  FUNCTION  6 

n 

In  this  appendix  we  shall  derive  and  briefly  discuss  an  expression 
for  the  form  of  the  function  0^  in  terms  of  other  dynamical  variables. 
To  simplify  the  notation,  explicit  mention  of  the  subscript  n  will  be 
suppressed  in  the  intermediate  steps.  The  starting  point  for  this 
derivation  is  with  the  second  and  third  of  equations  (5.10): 


A 

e. 


A 


p  =  sin  £  cos  0 


sin  £  sin  0  , 


whence  we  find  at  once  that 


tan  0(t) 


e2(t)  *  P(t) 

- K - 

ex(t)  •  p(t) 


(C-l) 


(C-2) 


Differentiating  this  expression  with  respect  to  time  and  rearranging 
terms  leads  to 


d0 

dt 


cos  0  r  />A  Ax  /A  \ 

"A  "A  2  t(erP)(e2'P> 
(e  L  *  p) 


A  A  A  A  A  Ax  A  As 

(e2*p)(e1*p)  +  (e1-p)(e2-p) 


/A  A\  /A  A\  1 
(e2*p)(e1-p)J 


(C-3) 
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2.  .  ,A  A  2  2 

But  cos  9/(e^*p)  =  cosec  Q  ,  and  the  use  of  several  vector  identities 

shows  that 


A  A  A  Ax  A  A  A  Ax  A  A 

(e1*P) (e2*P)  -  (e2-p)(eL.p)  =  (eL  x  e2) 


where  we  have  also  used  the  facts  that 
A  — ► 

k  •  oja  =  0  .  Therefore 
A 


A  -*  A 

p  =  U)  X  p  ,  and 


d0 


n  2  ,A  2  r  A  „  A  A  v  A  A  A  A 

=  cosec  £n(k.fi)  +  cosec  CnL  (e!  ‘  Pn>  (e2*Pn)  "  (e2*pn>  (el*  pn^  *  (C_4) 


This  division  of  9  into  two  terms  is  a  natural  one.  The  first 

n 

term  represents  the  contribution  to  5^  of  the  motion  of  the  target  and 
the  orientation  of  the  individual  dipole;  this  may  be  thought  of  as 
being  a  "dynamical"  term.  The  second  term  represents  the  contribution 
to  6n  of  the  change  in  the  coordinate  system  with  respect  to  which 

0^  is  measured,  due  to  the  motion  of  the  polarization  plane  itself  as 

A 

k  moves;  this  may  be  thought  of  as  being  a  "coordinate  transformation" 
term.  Unfortunately  the  second  term,  which  cannot  be  simplified  to  any 
significant  degree,  is  quadratic  in  one  of  the  unknowns:  p^  .  Therefore 
it  seems  very  probable  that  the  information  about  the  target  motion  and 
about  the  scatterer  orientation  which  is  contained  in  equation  (C-4) 
cannot  be  disentangled  and  extracted  from  the  equation.  Thus  the 
knowledge  of  0^  which  can  be  found  by  comparing  frequency  shifts  of 
the  diagonal  terms  of  the  transformed  scattering  matrix  will  have  to 
suffice  in  itself  and  not  as  a  prelude  to  more  useful  information  about 
the  target. 
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